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An updated heavy quark-antiquark potential is constructed by matching the short- 
distance perturbative part to long-distance lattice QCD results at an intermediate 
r scale. The static potential and the order 1/m potential are both analyzed in this 
, way. Effects of order 1/m in charmonium and bottomonium spectra are discussed 

in comparison. Charm and bottom quark masses are deduced from the spectra and 
C3 ' related to the quark masses of other schemes. 



I. INTRODUCTION 



>- , Investigations of the potential between a heavy quark and its antiquark have a long tra- 

dition. Early models featuring the basic Coulomb-plus-linear r dependence plus a hyperfine 
■ interaction established a quite successful phenomenology of charmonium spectroscopy below 
the open charm threshold. At short distances perturbative QCD is supposed to work. At 
larger distances lattice QCD results have confirmed the linearly rising confinement potential 
and established reliable values of the string tension. 

It has long been realized that the static quarkonium potential, constructed within the 
framework of QCD perturbation theory, has a badly convergent or even divergent behavior 
at short distances. This has been understood in terms of renormalon ambiguities. It has in 
fact been found that leading renormalon effects cancel in the sum of the static potential and 
twice the quark (pole) mass [1, 2], an important feature to be recalled later when we address 
the issue of quark masses in the context of the quarkonium potential and spectroscopy. 

The steadily increasing precision of current lattice QCD computations permits now an 
accurate matching of perturbative and non-perturbative approaches at an intermediate dis- 
tance scale. An update of this procedure, combined with an accurate analysis of charm and 
bottom quark masses, is the main topic of the present work. 

This paper is organized as follows. After a brief summary of first attempts in Section II, 
the perturbative static potential in coordinate space and its overall additive constant are 
derived in Section III and matched at an appropriate distance scale with results from lattice 
QCD. Section IV follows with an analogous construction for the potential at order 1/m in 
the quark mass. Section V discusses 1/m effects on bottomonium and charmonium spectra. 
In the last section charm and bottom quark masses are deduced and compared with existing 
determinations in other schemes. 



2 



II. PREPARATIONS 



As a starting point, consider the static quark-antiquark potential in momentum space at 
three-loop order: 

3 / ,,2 



+ y C(3)J CUT^n, 



where q is the three-momentum transfer. The coefficients a\ and a2 have been determined 
analytically and they read in the MS scheme [3-5]: 

31 20 

ai = -jC A - yT F n /? (2) 

/4343 A 2 7T 4 22 //o .\ 2 /1798 56 
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y - 16C(3)) C F T F n f + {jT F n^j , (3) 

where C F — 4/3, = 3, Tp = 1/2 for SU(3) and n/ is the number of light quark flavors. At 
three-loop order, infrared singular contributions ln(/z 2 R /<f 2 ) start to play a role (see e.g. [6]). 
The accompanying constant 

a 3 = 4 3 (209.884(1) - 51.4048n/ + 2.9061nJ - 0.0214nJ) (4) 

has been calculated independently in [7] and [8]. In this paper we focus on the two-loop 
level. In order to transform the potential to coordinate space, a s (\q |) in Eq. (1) is usually 
expressed (see e.g. [9]) as a powers series expansion in a s at some fixed scale //: 



a s (q) = a s (fx) 
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with t = ln(g 2 //i 2 ). The coefficients f3 n of the QCD (3 function are known up to four-loop 
order [10]: 

§^> = = _g, _ J| ft _ J^ft - + 0(o{). (6) 
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FIG. 1. Power series expansion of a s (q) at \x = 1 GeV for nj = 4. The thick line shows the full four- 
loop running according to the renormalization group equation. Note that the sign of a s {q — > oo) 
changes order by order. 



In this approach a Fourier transform leads to the standard, //-dependent definition of the 
coordinate space static potential, 



+ 



4tt 



3r 



+ 



47T 
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a 2 + p 2 (4gl(r) + n 2 /3) + 2< 7 »(2a 1 /3 +/3i) 



a 3 + 167r 2 Ci(ln(/i IH r)+7 E ) + /3 3 (8^(r) + 27rV(^) + 16C(3)) 



+ /3o(l2^(r) + 7r 2 )(a 1 /3 + 5/6/3 1 )+2^(r)(3a 2 /3 +2a 1 /3 1 +/3 2 ) + 0(a£) L (7) 



where g^r) = ln(/xr) + 7 E . The derivation of this r-space potential uses, in principle, 
information about a s (|<f|) over the full range in q space. However, the expansion (5) in 
powers of lng 2 is a good approximation only in a small neighborhood of the scale /i, as 
demonstrated in Fig. 1 for the choice ji — 2 GeV. Clearly, the behavior of a s (q) for q > 
10 GeV and q < 1 GeV is out of control for such an expansion. 

Fixing /i for instance at 2 GeV, the resulting coordinate space potential (7) behaves 
pathologically at r — > 0, as shown in Fig. 2. This behavior can be traced to the order- by- 
order sign changes observed in Fig. 1 for q— >oo. It can be improved using the renormalon 
subtracted scheme (see e.g. [11, 12]). Figure 3 shows the potential resulting from the fre- 
quently used ad hoc identification \i — 1/r which evidently works only at extremely short 
distances, r < 0.02 fm. Such a construction therefore does not suggest itself for a matching 
of the potential to lattice QCD results at typical distances r > 0.1 fm. 
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FIG. 3. Static r-space potential according to Eq. (7). Progressive orders are shown when fi is 
identified with 1/r. The choice ^ R = q 2 has been adopted at NNNLO. 



III. THE STATIC POTENTIAL 



Here we pursue a different strategy for constructing the static potential in coordinate 
space, based on the potential-subtracted (PS) scheme proposed by Beneke [1]. The r-space 
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FIG. 4. Static QCD potential (with rif = 3) from the restricted numerical Fourier transform (8). 
Shown is the NNLO potential for different values of (if. The curves have been shifted by a constant 
to match at small r values. 

potential is denned through a restricted Fourier transform as 

V ( °W/) = J j^e^V(°\\q\), (8) 
l<fl>M/ 

where V m (W\) 

is given in Eq. (1), but now a s (|<f |) for \q\ > \Xf is used without resorting to 
a power series expansion. The momentum space cutoff /if is introduced in order to delineate 
the uncontrolled low-g region from the high-g range where perturbation theory is considered 
to be reliable. The potential (8) differs from the "true" static potential, 

^(r)=/»(^) + V(°'(r 1 /i / ) ) (9) 

approximately by a constant, 

Pi 

v W M = J ^^V®m = ^fdq?Vt%) + Otf 1 *)^ (10) 
W\<Pf o 

which encodes non-perturbative low-g behavior that can be absorbed in the definition of 
the potential-subtracted (PS) quark mass (see Section VI). The correction of order ^"jr 2 is 
negligibly small in the range of interest (r < 0.1 fm). 

The potential V^(r, (if) is evaluated numerically using the four-loop renormalization 
group running of the strong coupling a s , see Eq. (6). For distances r < 0.2 fm, the resulting 
potential depends only marginally on (if as shown in Fig. 4. At the matching radius, 
r = 0.14 fm, the spread of (V^ — const) when varying (if between 0.7 GeV and 1.5 GeV is 
0.05 GeV. The convergence behavior of the potential is displayed in Fig. 5. Different orders 
have been matched at r = 0.01 fm and are then evolved to larger distances. Evidently, the 
convergence behavior of the potential V^(r, (if) is satisfactory. 
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FIG. 5. Static QCD potential (with rif = 3 and fj,f = 1.0 GeV) from the restricted numerical 
Fourier transform (8). Different orders of have been matched at 0.01 fm. The choice (J-i R = q 2 has 
been adopted at NNNLO. 
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FIG. 6. Static QCD potential for nj = 3, based on Eqs. (8) and (1), matched at intermediate 
distances to a potential from lattice QCD [15]. Dashed-dotted curve: simplest extrapolation using 
Coulomb-plus-linear r dependence. 

For bottomonium (n/ = 4 massless flavors), the input value for the renormalization group 
running of the strong coupling constant is chosen as a s (4.2 GeV) = 0.226 ± 0.003. In the 
case of charmonium (n/ = 3) we use a s (1.25 GeV) = 0.406 ±0.010 as input in the potential. 
These values are obtained from a s (mz — 91.1876 GeV) = 0.1184 ± 0.0007 [13] (for a theory 
with rif — 5 active quark flavors) and run down to 4.2 GeV and 1.25 GeV, taking into account 
flavor thresholds [14]. 
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The perturbative potential (8), valid at small distances, can be matched at intermediate 
distances to results from lattice QCD (see Fig. 6). We use a fit obtained by Bali et al. from 
a full QCD simulation [15]. The matching point (dashed line) is chosen at r = 0.14 fm. 
The exact position of the matching point is not important for the resulting shape of the 
potential. At r = 0.14 fm, both the perturbative and the lattice potential are expected 
to be reliable. Requiring that the first derivative of the potential is continuous at the 
matching point, we find for the cutoff in Eq. (8): /// = 0.908 GeV (bottomonium case) and 
fif = 0.930 GeV (charmonium case). The grey band reflects uncertainties in the Sommer 
scale r = 0.50 ± 0.03 fm (lattice part) and uncertainties in a s (\q |) (perturbative part) as 
given in the previous paragraph. This leads to a cutoff window: /i/ = 0.9^2 GeV (for both 
bottomonium and charmonium). The dashed-dotted line in Fig. 6 results from a simple 
Coulomb-plus-linear extrapolation from the lattice QCD data to short distances. Evidently, 
our more sophisticated perturbative QCD extrapolation based on Eqs. (8) and (1) differs 
from that simple form. 

For zero flavors one can check against accurate (quenched) lattice results [16] (see Fig. 7). 
Since a s cannot be extracted from experiment for rif = 0, we fit to the lattice points below 
0.12 fm. With a low momentum cutoff fif in the range 0.7-1.5 GeV, we find a s (1.25 GeV) = 
0.29 ± 0.01 for the flavorless strong coupling at the scale of the c-quark mass. The lattice 
scale r = 0.5 fm has been used here. A recent precision study of the zero-flavor case in a 
different approach can be found in Ref. [17]. 



IV. ORDER 1/M POTENTIAL 



The heavy quark- ant iquark potential can be expanded in inverse powers of the heavy 
quark mass m: 

K W=l A°> (r)+ ^ + yr) +0(1/m 3). 
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FIG. 8. The order 1/m potential with nj = 3 from the restricted numerical Fourier 

transform (13), for different cutoffs fi'j. The curves have been shifted by a constant to match at 
small r values. 



The perturbative potential at order 1/m in momentum space reads: 

y Wm = _ Ca + 0(a s )} , (12) 

^ \Q I 

with Cp = 4/3 and Ca = 3. This form is not unique, see Ref. [18], but we stick to the 
same convention as the one used on the lattice. V^^Qffl) can be transformed to r space as 
in Eq. (8), with a low momentum cutoff /x^ that may differ from \if. 

^ (1 W/)= / ^^V^(\q\). (13) 

\q\>n' f 

Evidently, the dependence of on the cutoff scale fi'j is again very weak for distances 
r < 0.2 fm as shown in Fig. 8. The variation of (V^ — const) when varying jx'j between 
0.7 GeV and 1.5 GeV is within 0.02 GeV 2 at the matching radius, r = 0.14 fm. This 
potential is again matched to corresponding results from lattice QCD [19, 20]. In order to 
fit the lattice data we have used 

V£ ] [r) = ~ + d! In (L^j + const , (14) 

with d = 0.0027 GeV 2 fm 2 , d! = 0.075 GeV 2 and an arbitrary length scale r$ that can 
be absorbed in the overall constant. The logarithmic form in Eq. (14) is motivated by 
effective string theory [21]. This parametrization extrapolates the lattice points for 
better than just a 1/r 2 + linear form. The lattice calculation of is quenched and subject 
to renormalization issues. A 15% uncertainty is therefore assumed in the lattice potential, 
in addition to the uncertainties in the Sommer scale tq. At short distances a deviation of 
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FIG. 9. The order 1/m potential V^(r, fi'*) with rif = 3, based on Eq. (13), matched at in- 
termediate distances to a potential from lattice QCD [20]. The fit and the matching has been 
performed using the lattice QCD points with /3 = 5.85 (circles). Legends are analogous to Fig. 6. 
Dashed-dotted curve: simple extrapolation using Eq. (14). 



the perturbative potential V^(r, n'j) from \Q t of Eq. (14) (dashed-dotted line in Fig. 9) 
is apparent. For the cutoff in Eq. (13) with error estimate we find: fi'^ 
(bottomonium case) and fi'j = 1.6~^ 8 GeV (charmonium case). 
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V. QUARKONIUM SPECTRA 

Given the potential V = V {0) + 2 V {1) /m up to order 1/m in the heavy quark mass, we 
can now examine the resulting bottomonium and charmonium spectra, with focus on the 
effects of the 1/m term. The Schrodinger equation 

_ ^!v 2 + 2m rs (/i /)/ / / ) + vM{r,n f ) + -V^(r^' f ) - E 
m 1 m 1 

is solved with the fixed values for fif and jjilf as derived in the construction of the poten- 
tials. The potential-subtracted mass mg§(///,^), to be defined and discussed in detail in 
Section VI, is the only free parameter. It sets the overall energy scale and it is ultimately 
fixed by comparison with the measured bottomonium and charmonium spectrum. 

The (heavy) quark mass is not directly measurable. The mass m appearing in the de- 
nominator of the kinetic energy and the term is therefore not a priori determined. 
In practice we use values close to the static MS masses: m = 4.2 GeV for bottomonium 
and m = 1.25 GeV for charmonium. Small variations from these values do not have any 
significant influence on the results. 

Consider now first the bottomonium spectrum below BB threshold 1 (see Fig. 10). One 
can start by fixing mgg(///,/ij) such that the measured T(2S) energy is reproduced. Alter- 

1 Above the BB threshold the bb potential develops an imaginary part and the present strategy (including 
lattice QCD) does not apply. 



^(f) = 0, (15) 
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FIG. 10. Bottomonium spectrum in comparison with experiment. Static plus order 1/m results are 
shown on the left, with additional hyperfine effects (h.f.) added phenomenologically using Eq. (16). 



natively, the center of the Xb(lP) triplet can be used for calibration. These states remain 
almost unchanged by the 1/m effects. The more tightly bound r/b(lS) and T(1S) states 
respond, as expected, more sensitively to the non-static corrections induced by V' 1 '^) with 
its pronounced behavior at short distances. 

An additional effective one-gluon exchange spin dependent term, 



+ 



at s ( L ■ B\ + L ■ <72 



, (16) 



with af = 0.3 would move all IS and IP states well into their observed positions. For this 



11 



> 
o 





3.8 


> 


3.7 






o 


3.6 


CO 




CO 


3.5 


s 


3.4 




3.3 




3.2 



D_D-threshold_ 

77c (2S) 




JAK1S) 



static 


+l/m 


+h.f. 


experiment 






DD-threshold 


: IP 




h c (lP) Xcj(1P) : 
















i Pl 3 Pj ; 


static 


+l/m 


+h.f. 


experiment 



FIG. 11. Charmonium spectrum in comparison with experiment. Static plus order 1/m results are 
shown on the left, with additional hyperfme effects (h.f.) added phenomenologically. 



purpose we replace the delta function (that is exclusively sensitive to the wave function at 
r = 0) by a Gaussian distribution 



6®( 



r -> 



-r 2 /a 2 



(17) 



with a = 0.03 fm. Of course, these procedures based on Eq. (16) are purely ad hoc and 
need to be substituted by the full potential of order 1/m 2 , to be investigated in forthcoming 
work. 

As expected, the influence of the 1/m term in the potential is much stronger for char- 
monium than for bottomonium (see Fig. 11). For charmonium we choose the potential- 
subtracted mass in Eq. (15) (reflecting the unknown constants in and V^) such that 
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the measured h c (lP) energy is reproduced without 1/m corrections. With this choice, how- 
ever, the part produces a downward shift of 127 MeV in the IS states (r/ c and J/0) 
relative to the static result. This shift is too large in comparison with the measured r/ c and 
J/0 energies. 

The IS states are naturally more sensitive to 1/m corrections than IP and 2S states 
because of the leading 1/r 2 short-distance behavior of V^. Hence the large shift of the 
IS energy level at order 1/m does not come unexpectedly. It is nevertheless evident that, 
no matter which choice is adopted for adjusting the unknown constant in the IS and 
IP states of charmonium cannot be simultaneously reproduced at order 1/m. Unlike the 
situation in bottomonium, corrections of order 1/m 2 are presumably large in charmonium. 

A manifestation of substantial 1/m 2 effects is the relatively large observed splitting of 
117 MeV between r/ c and J/0, driven by an effective coupling strength (a| ff /m 2 with af 1 = 
0.3 in the phenomenological SV spin of Eq. (16)) that is an order of magnitude larger than 
for bottomonium. A systematic investigation of the 1/m 2 potential, in Eq. (11), is 
mandatory now that lattice QCD data for are becoming available [22, 23]. 



VI. CHARM AND BOTTOM QUARK MASSES: 
POTENTIAL-SUBTRACTED AND MS SCHEMES 



The static potential is determined up to an overall constant. Introducing the potential- 
subtracted (PS) quark mass as 

d 3 q 



mps(pf) = m polc + - 



1 4^ mm=m ^ + ^ S iqq2Vm(q) (i8) 

\q\<Hf 

in terms of the pole mass of the charm or bottom quark, this mps absorbs the unknown 
constant and does not suffer from the leading renormalon ambiguity [1, 2]. With the previous 
choice of the constant to reproduce the measured T(2S) energy, we find the value mps(/// = 
0.908 GeV) = 4.78 GeV in the PS mass scheme. In the charmonium case we fit to the h c (lP) 
energy and find m PS (/i/ = 0.930 GeV) = 1.39 GeV. To convert the PS mass to the mass in 
the more commonly used MS scheme, it is necessary to introduce the pole mass m po i e as an 
intermediate step. The relation between m po i e and mps reads [1]: 
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"Vic = mps(ttf) + 



C F a a (n)nf 
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«2-(2«j >„- ; In ^ 
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2) f ^Mn 2 ^-41n^ + 



fi 2 



+ 0(a 3 s )\, (19) 



with the same conventions as in Section II. Note that a renormalization scale \x appears in 



the coupling. In the following fj, is set equal to the MS mass m 



m Ms( m MS ) • This m is not 
known at that point and has to be computed iteratively. In a second step the pole mass is 
converted to the MS mass [24, 25]: 
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FIG. 12. Scale dependence of mps([if) for the bottomonium case. The solid curve shows the 
numerical [if dependence from Eq. (18). The dashed and dashed-dotted curves show the /in- 
dependence given by Eq. (19) at different matching points. Point (2) is preferred to point (1) for 
the determination of my,. 

Note that both relations (19) and (20), taken individually, show a poorly convergent behavior 
whereas the relation between mps(fif) and m is expected to be stable. This is in fact 
confirmed numerically. Using [if = 0.908 GeV (from Section III) one obtains the value 
% = 4.27 GeV for the bottom quark; with [if — 0.930 GeV one finds m c = 1.24 GeV for 
the charm quark. 

However, as shown in Fig. 12, the [if dependence of mps([if) (Eq. (19)) differs for /if <Cm 
from the [if dependence coming from variations of the cutoff in the numerical integral (solid 
line). Instead of matching at [if = 0.908 GeV it is obviously preferable to use the numerical 
[if dependence first to translate mps (/•«/) into mps(m) and then apply Eqs. (19) and (20) 
to translate this value into the MS scheme. By this method the extraction of m becomes 
independent of the value of [if used in the construction of the potential. This leads to 
improved mass values, Mb = 4.20 GeV for the bottom quark and m c = 1.23 GeV for the 
charm quark in the MS scheme. 

To determine values for the quark masses at order 1/m a redefinition of the PS mass is 
required: 

m^{[if,[i'f) = mpsO/) - —C F C A a 2 s (fn) [if . (21) 

The 1/m-term stems from an analogous calculation as in the static case. The renormalization 
scale [i that appears in the coupling has again been identified with m. We determine 
mfs([tf, [t'f) for the b- and c-quark by fitting to the empirical T(2S) and h c (lP) energies, 
respectively, and convert these values numerically to m^im, m). This leads in a second step 
to the 1/m-improved MS values m b = 4.18 GeV for the bottom quark and m c = 1.28 GeV 
for the charm quark. In Table I the quark masses in the MS scheme found in our approach 
are summarized and compared to values given by the Particle Data Group [26]. We have 
performed error estimates for the quark masses, reflecting uncertainties in the potentials 
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MS masses [GeV] 




Static 


Static + 0(l/m) 


PDG 2010 


Bottom quark 
Charm quark 


4.20 ± 0.04 
1.23 ±0.04 


A 1S+0-05 
i no+0.07 


4 iq+0.18 

1 27 +ao7 

L - Zl -0.09 



TABLE I. Comparison of quark masses obtained in our approach (leading order plus order 1/m 
corrections) with the values listed by the Particle Data Group (PDG) [26]. See text for details 
concerning error estimates. 

(static and order 1/m). Additional uncertainties are included from our specific choice of 
matching to the empirical T(2S) and h c (lP) energies. The errors at order 1/m have increased 
in comparison to those for the static case since they incorporate in addition the error band 
from The error estimates at order 1/m do not include possible further uncertainties 

appearing at order 1/m 2 . 

VII. SUMMARY 

Improved bottomonium and charmonium potentials have been derived up to and includ- 
ing order 1/m in the heavy quark masses by systematically matching perturbative QCD 
results to accurate lattice QCD data at an intermediate distance scale, r = 0.14 fm. A sin- 
gle constant (the potential-subtracted quark mass) is adjusted to reproduce the T(2S) and 
h c (lP) masses, respectively. The predicted pattern of all other subthreshold states agrees 
well with the empirical bottomonium spectroscopy. For charmonium, 1/m effects are far 
more pronounced, as expected. The potential-subtracted heavy quark masses at order 1/m, 
when translated to MS masses, agree well with those listed by the Particle Data Group. In 
a next step, corrections of order 1/m 2 will be studied along the same lines. 
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